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An adaptive differential pulse-code-modulation technique for en- 
coding and decoding has been implemented using the Bell Labora- 
tories digital signal processor integrated circuit. The encoder/de- 
coder operates in real time and can accommodate 3- or 4- bit (8 kHz) 
encoding. In this paper, we discuss details of the implementation, the 
basic algorithm, and the features utilized in the digital signal proc- 
essor. 

I. INTRODUCTION 

Adaptive differential pulse-code-modulation (adpcm) encoding has 
been shown to be a simple and effective method for digitally encoding 
speech at bit rates in the range of approximately 24 to 48 kb/s. 1 " 6 At a 
rate of 24 kb/s, adpcm can provide a good quality reproduction of 
speech that is acceptable for applications such as computer-controlled 
digital voice response systems. 4 At a rate of 32 kb/s, it can provide 
essentially a telephone bandwidth "transparent quality" (a quality 
that is indistinguishable from the original uncoded source) for a single 
tandem encoding. Adaptive differential pulse-code modulation has 
been studied for use in some types of transmission systems, 6 and for 
message storage and retrieval systems in which a reduction by a factor 
of two in bit rate over that of conventional 64 kb/s /i-law companded 
pcm is desired. 4 

Various forms of hardware have been suggested for the implemen- 
tation of adpcm coders. Some designs are based primarily on analog 
hardware 1 ' 3 where parameters are pair- wise tuned between transmit- 

1547 



ters and receivers. This results in problems with repeatability and 
stability in the analog designs. Bates 5 and Adelman, Ching, and Gotz b 
subsequently presented two different methods of designing all-digital 
adpcm coders. The Bates approach uses a ttl logic design with a ROM- 
based look-up table for the adaptive step-size and an up/down counting 
scheme with digital adders and subtracters to avoid the use of a digital 
multiplier. The hardware requires a twelve-bit linear pcm input and 
has a total "package count" of approximately 80 standard ttl logic 
packages. It was constructed on two wire-wrapped Augat cards (one 
encoder and one decoder per board). The hardware by Adelman et al., 
was a ROM-based design that accepted a standard 64-kb/s ,u-law com- 
panded pcm signal. 

Recently, an LSI digital signal processor (dsp) has been developed 
by Bell Laboratories. 7 The DSP is a programmable processor capable 
of performing the entire adpcm algorithm for multiple channels in a 
single lsi device. The ability of the processor to convert between 
conventional /i-255 companded pcm and two's complement binary code 
formats allows the adpcm algorithm to use either format. The config- 
uration of the adpcm algorithm that is implemented on the dsp is 
described in Section II. Section III discusses details of how the algo- 
rithm is configured to use dsp features. Sections IV and V describe the 
hardware configuration and measured performance, respectively. 

II. THE ADPCM ALGORITHM 
2. 1 Overall configuration 

Figure 1 illustrates the basic configuration of the dsp implementa- 
tion of adpcm. The input analog signal s(t) is sampled and a/d 
converted to an 8-bit /u-law companded pcm format to produce the 
sampled data signal s(n), where n is the discrete time index. These 
operations are done externally to the dsp. Then s(n) is converted in 
the dsp from /i-law format to a 20-bit linear pcm format for internal 
processing. 

The adpcm encoding is performed entirely within the transmitter 
dsp. The output is a 3- or 4-bit codeword, I(n), which can be obtained 
through the normal output channel of the dsp. In the receiver, a 
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Fig. 1— Overall configuration of the adpcm encoder/decoder using the dsp. 
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second dsp is used to decode I(n) back to an 8-bit /x-law format. This 
signal is then converted to an analog signal with an external d/a 
converter. The clock and framing signals are carried between trans- 
mitter and receiver on separate lines. The transmission protocol for a 
given application depends on the data channel in which the system is 
implemented. 

2.2 Block diagram of ADPCM 

Figure 2 shows the block diagram of the adpcm algorithm that was 
implemented on the dsp. The transmitter is shown in Fig. 2a and the 
receiver is shown in Fig. 2b. We assume that s(n) in Fig. 2 is in a 20-bit 
linear pcm format, because of a direct 20-bit input to the dsp or 
because of an 8-bit /i-law to 20-bit linear pcm conversion that has been 
performed within the dsp. 

The adpcm algorithm can be partitioned into three basic parts (see 
Fig. 2): the adaptive pcm quantizer, the differential predictor loop, and 
the adaptive gain (step-size) control for the quantizer. We discuss 
these operations in the above order since this is the order in which 
they are computed. 

2.3 Adaptive quantizer 

A predicted value, p(n), is first subtracted from the input signal s(n) 
to form the difference signal 

e(n)=s(n)-p(n). (1) 

The value of p(n) is obtained from the predictor and is based on 
computations performed at the previous sample time, n — 1. 

The difference signal e(n) is then quantized by the adaptive quan- 
tizer to produce the adpcm codeword I(n). This adaptive quantization 
is achieved by first scaling e(n) to the range of a 3- or 4-bit fixed step- 
size quantizer (see Fig. 2a). The scaled signal is denoted as 

eAn) = V(n).e(n), (2) 

where V(n) is an adaptive scale factor that is also determined from 
data available at the previous sample time n — \. The combination of 
scaling by V(n) followed by the fixed step-size quantizer is equivalent 
to a quantizer with an adaptively varying step-size (which is inversely 
proportional to V(n)). 

Figure 3 shows the characteristics for a 3-bit (8 level) fixed step-size 
quantizer. The input signal e a (n) (consisting of an integer plus a 
fractional part) is converted to one of eight quantization levels in the 
range —7/2 to +7/2 corresponding respectively to integer codewords, 
/(n), in the range —4 to +3. The output quantized signal is denoted as 
e a (n) and is related to the code word I(n) by the relation 
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Fig. 2b — Block diagram of the adpcm decoding algorithm. 



e s (n) = I(n) + 0.5. 



(3) 



A similar quantizer characteristic is used for the 4-bit (16 level) design. 
The unsealed, decoded difference signal can be obtained from e s (n) 
by rescaling it by A (n), where A(n) is inversely related to V{n) and is 
directly proportional to the "step-size" of the equivalent adaptive 
quantizer. Thus, 



e(n) = b(n)-e s (n) 



(4) 
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Fig. 3 — Quantizer characteristic for the 3-bit pcm quantizer. 
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is the adaptively quantized representation of the difference signal e(n). 
This completes the adaptive quantization part of the algorithm. 

2.4 First-order predictor 

The sum of p(n) and e(n) also forms the adaptively quantized 
representation of the input signal (see Fig. 2a) 

s(n) = p{n) + e(n). (5) 

However, the predictor value p(n + 1) and the quantizer scale factors 
V(n + 1) and A(/i + 1) need to be updated for the next sample time 
n + 1. 

The new predictor value is computed as fi times s(n); i.e., 

p(/i+l)=y8.s(n), (6) 

where /? is the predictor "leak" factor. 2 A value of (S = 0.85 is often 
used for speech encoding. The value of p(n + 1) is then stored for use 
at the next sample time n + I. 

2.5 Adaptive step-size control 

The computation of the new value of A(n. + 1) and V(n + 1) for the 
next time sample n + 1 is more involved. The algorithm used here is 
based on the robust step-size adaptation approach. The details and 
advantages of this form of the algorithm are discussed in considerable 
detail in Refs. 8 to 10. The method is as follows. The step-size referred 
to above is the effective spacing in the quantizer levels observed in the 
unsealed quantized signal e(n). Since this spacing is proportional to 
A(rc), we will refer to A(w) in the following discussion as the step-size. 

In the robust adaptation algorithm, the new step-size A(/i + 1) is 
chosen as 

A(n + l) = (A(n)y.M(I(n)). (7) 

That is, it is the old step-size A(/i) raised to the power y (0 < y < 1, 
typically y = 0.98) and scaled by a factor M( ■ ) which is a function of 
the code word I(n). If the code word is one of the upper magnitude 
levels of the quantizer [e.g., I(n) = -4, -3, 2, or 3 in Fig. 3], a value of 
M( ■ ) greater than one is used to increase the step-size of the quantizer 
for the next sample time. If the codeword is one of the lower magnitude 
levels [e.g., I(n) = -2, -1, 0, 1 in Fig. 3], a value of M{- ) less than one 
is used to reduce the quantizer step-size for the next sample time. In 
this way, the algorithm continually attempts to adapt the step-size 
such that the dynamic range of the adaptive quantizer is matched to 
the range of the signal e(n). 

A more direct approach for implementing this algorithm is obtained 
by expressing eq. (7) in logarithmic form. Let 
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rf(#i) = log(A(n)), (8) 

and 

m(I(n)) = log(M(I(n))), (9) 

where the base of the logarithm is determined by the choice of 
parameters in the coder, as discussed later. Then, eq. (7) becomes 

d(n + 1) = yd(n) +m(I(n)), (10) 

and it has the form of a first-order difference equation. This is the 
equation that is implemented by the upper right-hand part of the 
block diagram in Fig. 2a. 

The driving function m(I(n)) is a function of the code word I(n), 
and it is determined in the adaptation logic algorithm which performs 
the function 

IT( u / 8 if \I{n) +0.51 >2.5, ,„„ , 

m(/(/l)) = i-3if|/(n) + 0.5|<2.5 (lla) 

for the 3-bit quantizer and 

ITt v . f 8 if \I(n) + 0.51 > 4.5, 
m(/(n)) = |_ 3if j /(/i)+05 ] <45 ' (lib) 

for the 4-bit quantizer. 

The output of the step-size adaptation loop is the signal d(n), which 
eq. (8) shows as the logarithm of the desired step-size. 

Thus, to obtain A(/i) (and V(n)), we need to implement the relations 

A(n) = exp(d(n)) (12a) 

and 

V(/l) = AM = ex P (-rf(,l)) ' < 12b ) 

where, again, the base of the logarithms and exponentials in eqs. (8) to 
(12) are determined by the choice of parameters of the coder, as 
discussed later. 

The exponentials in eqs. (12a) and (12b) are computed by using 
look-up tables that are stored in the program rom as indicated in Fig. 
2a. The integer value of d(n) is taken and constrained to the range 
-32 < [d(n)] < 32 for a lookup table size of 64. Table offset values are 
then added to this value to produce the appropriate rom address 
locations. The tables are set up so that a value of d(n) = points to 
the center of both tables (V[32] and A[32]). Note that the bracketed 
values V [•] and A [•] in the tables of Fig. 2a refer to the contents of 
table locations addressed by [•] and not to sample times (which are 
referred to by the parenthesis notation V(rc) and A(/i)). 

The above algorithm uses 64 different step-size values of A(/i) and 
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V{n). The table values are chosen to span the desired dynamic range 
of the input signal. If signals of larger or smaller amplitude than this 
are encountered, the step-size saturates at the upper or lower table 
values, respectively. 

The table address locations are stored for use at the next sample 
time n + 1 to access values A(n + 1) and V(n + 1) as they are needed. 
This completes the operation of the adpcm encoder. 

2.6 Parameter selection for the step-size tables 

The parameters of the step-size table were chosen such that the 
ratio of the maximum to minimum step-sizes in the table is 256, i.e., 

A[64]/A[l] = 256 = 2 8 

and 

V[l]/V[64] = 256. 

This gives a step-size adaptation range of 48 dB (8 bits) that is 
appropriate for speech coding. Since there are 64 logarithmically 
spaced step-size values in each table, this corresponds to a 0.75-dB 
resolution, i.e., step-sizes increase by a ratio of 1.0902 in the table, 

A[i] = A[l]. (1.0902) i ~\ 

V[i] - V[l]-( 1.0902) _, ' +1 . 

The maximum step-size is chosen so that the maximum range of 
e(n) (approximately 17 bits) scales to the maximum quantizer range 
(3 or 4 bits). This prevents overloading on the high end of the dynamic 
range. Thus, for a 4-bit quantizer 

A(/i)| max = A[64] = 2 17 - 4 = 2 13 

and 

A(«)| mi n = A[l] = 2 5 . 

For a 3-bit quantizer, values of A[ • ] should be increased by two and 
values of V[ ■ ] should be decreased by two for best dynamic range 
performance. 

The manner in which these table values are scaled and stored in the 
dsp will be explained more fully in Section 3.4. 

2.7 Decoder for ADPCM 

Figure 2b shows a similar block diagram for the adpcm decoder. The 
input code word I{n) is converted to the decoded difference signal e(n) 
by adding 0.5 (see Fig. 3) and scaling the result by A(n). The decoded 
output signal s(n) is then obtained by accumulating values e(n) in the 
predictor loop in the same manner as in the encoder. The new step- 
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size is then computed in an exactly duplicate manner to that in the 
encoder. Thus, the adpcm decoder duplicates a subset of the encoder 
block diagram. 

III. PROGRAMMING TECHNIQUES UTILIZED IN THE DSP FOR THE 
ADPCM ALGORITHM 

The adpcm algorithm, as it is configured above, is designed so that 
it can be conveniently implemented on the dsp. We have already 
alluded to ways in which the dsp is used in implementing the algorithm. 
In this section, we point out some additional aspects of the program 
and discuss how some of the unique features of the dsp are used. 
Discussions in this section use the 4-bit algorithm as an example. 

3. 1 Memory utilization 

The encoder program and the step-size tables for A[ • ] and V[ • ] are 
stored in rom occupying approximately 170 words of memory. Five 
ram locations are used for storing the following values: 2(1 {n) + 0.5), 
2 s(n), p(n), d(n) and the integer version of d(n). Access to the step- 
size tables is made by setting the rx register in the dsp to the center 
address of the desired table. The table address for the appropriate 
step-size is then obtained by setting the k register to the integer value 
[d(n)] (limited to the range —32 to 31) and then incrementing rx by 
the value k using the rxk command. The rx register then contains the 
rom address for the appropriate step-size V[.]. The dsp instructions 
that implement this technique are as follows: 

rya = 5; "ram pointer to int [d(n)]" 

rx = &TABLE; "set rx pointer to V[32]" 

k = rym;;; "set k to int [d(n)]" 

a = pp = rxk*c; "pointer to appropriate V[-]" 

Note that this technique only works for table sizes up to 256, since the 
k register is limited to 8 bits. 

3.2 Quantization 

The computation of the quantizer input, e s (n), is accomplished in a 
straightforward manner on the dsp by using a subtraction and a 
multiplication. The conversion from e s (n) to a 4-bit integer I(n) 
(according to the quantizer characteristic in Fig. 3) is done by truncat- 
ing the fractional part of e s (n). 

Assuming e(n) has been computed and stored in the w register and 
rx is pointing to the appropriate step size, the following dsp instruc- 
tions compute I'(n). 
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p = 8*c; 
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if (a<0) doau (); 
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p = l*c; 

a = p p = rxj*w; "compute e s (n) = V(n)-e(n)" 

a = p + a/2 p = — l*c; "compute e s (n) + 0.5" 

a = p & a; "truncate to form Fin)" 

The truncation is achieved by using the accumulator control statement 
a = p & a, which performs a bit-by-bit AND operation between the p 
and a registers in the dsp. T{n) is stored in the a register and the 
(binary) number 1111 • • • 111.000 • • • is stored in the P register by the 
instruction p = — l*c (c = 2 14 ). The operation zeros out the fractional 
part of (e s (n) + 0.5) and leaves the integer part untouched. The 
resulting integer I'(n) is then constrained to the range —8 < I'(n) < 7 
(for 4-bit quantization) to form the desired codeword I(n). This is done 
using the conditional au operation as shown in the following instruc- 
tions. 



"if (F(n) < -8) T(n) = -8" 

"if (/'(/i)> 7) /'(«) = 7" 
"4-bit I(n)" 
This yields the 4-bit value I(n) in the accumulator. 

3.3 Internal scaling of data 

At the output of the quantizer a value of 0.5 is added to I(n) (see 
Fig. 2) to produce e s (n), which is then multipled by A(n). To accom- 
plish this, e s (n) must be transferred from the a register to the w 
register and the input of the multiplier. Since this transfer from the a 
to w registers truncates the fractional part of the value in the a 
register, it is first scaled by a factor of 2 to avoid the truncation of the 
fractional part. This scale factor of two is carried through the compu- 
tation of s(n) in the block diagram of Fig. 2. 

3.4 Scaling of step-size table values 

The values stored in the step-size tables in rom are more conven- 
iently handled if they are scaled to be less than 2 in magnitude. The 
internal dsp arithmetic is such that 16-bit values from ROM are assumed 
to have 14 fractional bits. This is appropriate for the inverse step-size 
table V[-] in which numbers (e(n)) in a 17-bit range are scaled down 
to a 4-bit range (e s {n)). This table takes on values: 
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V[l] = 0.031250 
V[2] = 0.028617 

V[63] = 0.000133 
V[64] = 0.000122 

which can be directly stored in the ROM. 

In the step-size table A[.], the inverse of these numbers must be 
stored. To accomplish this, the inverses are first scaled by a factor of 
2 -14 and, thus, they take on the fractional range (denoted by primes): 

A'[l] = 0.001953 

A'[2] = 0.002133 

A'[63] = 0.457650 
A'[64] = 0.499754 

When these table values are used in the multiplication, the resulting 
product A'(/i)e s (/i) is scaled back up by a factor 2 14 using the 14-bit 
shift option a = a «: 14. 

The inverse relationship between the table values requires that 

V[i].A'[i].2 14 = l, 
*=1,2, ... ,64. (14) 

Since V[i] and A'[i] must be quantized to 16-bit numbers (14-bit 
fractions), for storage in the rom, the condition in eq. (14) cannot be 
met exactly. To obtain the greatest accuracy in representing these 
numbers, the smaller of the two values V[i] or A'[i] for each value of 
i is quantized first. The reciprocal of this number is then computed 
(with floating point accuracy) and scaled by 2 -14 . This value is then 
quantized to the 16-bit range of the rom to produce the inverse table 
value. Thus, the accuracy of the condition in eq. (14) is maintained as 
closely as possible. 

3.5 Control of the address range for the tables 

Another unique part of the program involves the control of the range 
of the table address pointer d(n) to the range -32 to 31 (excluding the 
table offset). This is accomplished with the aid of the overflow protec- 
tion feature of the dsp when data is transferred from the a register to 
the w register. The operations in the step-size control loop are scaled 
so that d(n) is computed in the upper range of the a register. When 
this number is transferred to the w register, it is automatically limited 
to the proper range by the overflow protection option in the dsp. 
Scaling this number back down to a 6-bit integer range gives the 
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desired range of -32 to 31 for the table location (excluding the table 

offset). 

Specifically, this is accomplished as follows. First the driving func- 
tion 

m w \ 8.2 8 if|2(/(/i) + 0.5)|>9, 
m(I{n)) - |_ 3 . 2 8 if | 2 (/(n) + 0.5) | < 9 

is computed and stored in the a register. The value m(I(n)) is then 
multiplied by 8 twice and added to yd(n) to form d{n + 1) scaled by 
2 14 . The value d(n + 1) is limited when it is transferred to the w 
register, and the table look-up offset, int (d(n)), is obtained by multi- 
plying the w register by 2 10 . This result is directly loaded into the k 
register as discussed in Section 3.1. Assuming the a register contains 
the absolute value of 2 (I(n) + 0.5) - 9, and the RY register points to 
d(n), the following dsp instructions compute d(n + 1) and the table 
offset value, integer d(n). 

p = 04000»c; 

if (a < 0) doau ( ); "if (a > 0) m = 8*2 8 " 

a = p p = 0176400*c; "else m = -3*2 8 " 

a = p p = 0»c; 

a = p + 8»a p = .98»rym; "compute yd(n)" 

a = p + 8»a p = 0*c; "compute d(n + 1)«2 



w = a; 



"overflow protection" 



rdp = w a = p p = 16*w; "save d(n + 1)" 

a = p; 
w = a; 
rdp = w; "save int(d(n))" 

IV. HARDWARE CONFIGURATION 

The hardware used to implement the algorithm described in this 
paper consists of two 16.5- by 11.5-cm wire wrap cards, one for the 
transmitter and one for the receiver. Both cards are of identical 
construction with jumper plugs used to determine whether the card 
will act as a transmitter or a receiver in the configuration shown in 

Fig. 1. 

Figure 4 shows a more detailed block diagram of the cards. Each 
card contains a dsp, along with a 2048-word by 16-bit prom holding 
the program and step-size tables. (For permanent applications the ROM 
would be integral to the dsp). In addition, there are analog filters, \i- 
law encoder and decoder chips, clock generators, and synchronization 

logic. 

The adpcm transmitter card accepts an analog input which is applied 
to a buffer amplifier and then to a low-pass filter and a /i-law encoder. 
The sampling rate of this encoder is determined by the repetition 
period of the sync signal produced by the sync generator. The serial 
output of the encoder is shifted out by the clock signal. These three 
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signals (data, sync, and clock) are connected through a jumper plug to 
the data in (di), input sync (isy), and input clock (ick) pins of the dsp. 
Operating from the encoding program stored in the prom, the dsp 
converts the incoming data into adpcm code words which become 
available at the data output pin (do). The shifting of this word out of 
the dsp is controlled by output sync (osy) and output clock (ock) 
lines. These three signals are connected through another jumper plug 
to the external output cable which sends them to the receiver card. 
The ju-law decoder on the transmitter card is not used. 

The adpcm receiver card receives the data, sync, and the clock 
signals from the transmitter card over the external input cable, which 
is connected through the first jumper plug directly into the dsp. The 
ju-law encoder and sync generator on this card are bypassed. The dsp 
takes the incoming adpcm code words and converts them to jt-law 
pcm, according to the instructions of the decoding algorithm in its 
prom. The jt-law data words are shifted out of the dsp on the do line 
and are connected through the other jumper plug to the input of the 
/i-law decoder chip. The analog output of the decoder is bandpass 
filtered and then sent through a buffer amplifier to produce the 
reconstructed analog output signal. 

The cards used here are not limited to only one application of the 
dsp. By setting the jumper plugs so that audio input and output are 
both on the same card, different types of filtering algorithms can be 
tested. Conversely, both the input and output of a card can be con- 
nected over the external cables. In this manner, several dsps can be 
connected together serially for more complex signal processing. 

V. PERFORMANCE 

Figure 5 shows the signal-to-noise ratio measured for the 4-bit 
ADPCM//i-law coder (Fig. 1) as a function of frequency for input signal 
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levels of 0, —20, and —40 dB of full scale. It corresponds closely to that 
of the coders in Ref. 3. The dynamic range is limited to about 48 dB 
because of the range of the step-size tables and the /x-law encoders. 
This dynamic range performance can be modified (within the limita- 
tions of the ju-law encoders) by changing the step-size tables. At low 
input levels the significant noise is that of the /i-law encoders. 

The subjective quality of the 4-bit adpcm is very similar to 200- to 
3200-Hz (telephone band) filtered speech without encoding. This sug- 
gests that a single adpcm link can provide essentially a "transparent" 
quality for telephone bandwidth speech. 

VI. CONCLUSIONS 

In this paper, we have discussed the implementation of the adpcm 
algorithm on the dsp. The encoder program uses 26 percent of the 8- 
kHz real-time capability of the dsp running with a 5-MHz clock. It 
uses 4 percent of the ram and 17 percent of the ROM. The decoder 
program uses 22 percent of the real-time capabilities, 4 percent of the 
ram and 10 percent of the rom. This suggests that 4 encoders, 4 
decoders, or 2 encoder-decoders could be implemented on a single dsp. 
The program uses a number of unique features of the dsp to achieve 
an efficient implementation of the algorithm, and it demonstrates the 
flexibility of the dsp in doing small-to-medium scale algorithms of this 
type. 
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